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Abstract 

Tha impact of spacecraft acceleration environment on 
Protein Crystal Growth (PCG) is studied. A brief overview 
of the Space Shuttle acceleration environment is provided 
followed by a simple scaling procedure used to obtain 
estimates of the flow and concentration field characteristics 
in PCG. A detailed two-dimensional numerical model is 
then used to simulate the PCG system response to different 
disturbance scenarios; viz. residual g effects, impulse type 
disturbances and oscillatory inputs. The results show that 
PCG is susceptible to g— jitter and is a good candidate for 
vibration isolation. 

Introduction 

Crystal growth in space benefits not only from its 
reduced gravity environment but also from the absence of 
the hydrostatic pressure which assists certain crystal growth 
and refinement methods. Gravity driven phenomena are 
thus reduced in strength, and a purely diffusive fluid 
behavior can be attained. I n addition, past materials science 
experiments have shown that microgravity can also help 
produce larger crystals. While gravity related effects are 
definitely curtailed in space, they are nevertheless present 
to some degree due to the acceleration environment on 
board the spacecraft These disturbances, collectively 
referred to as g-jitter, can be deleterious to certain 
experiments where the minimization of the acceleration 
level is important Advanced vibration isolation techniques 
can be utilized to actively fitter out some of the detrimental 
frequencies and help in obtaining optimum results. 
However, the successful application of this technology 
requires the detailed analysis of candidate fluids 
experiments to gauge their response to g-jitter and to 
determine their acceleration sensitivities. 

Several crystal growth experiments in the Protein 
Crystal Growth (PCG) area, besides others, are expected to 
be carried out on future shuttle flights and on the space 
station. The need for vibration isolation systems or 
components for microgravity science experiments can be 
to grow as experiments and available hardware 
become more complex This technology will also And 
increased application as the science community develops 
an awareness of their specific needs relative to the 
environment available In manned space missions. 
Vibration isolation research strives to develop a 
microgravity environment requirement that defines 
tolerance limits on the allowable g-level, and provides the 
required technology to achieve it This effort will assist in 


establishing the tolerable acceleration levels for specific 
experiments. This paper is directed towards modeling PCG 
with the following objectives in mind. 

I. To computationally determine vibration sensitivity of 
Protein Crystal Growth Experiments. 

2 Determine if these experiments can benefit from 
vibration isolation techniques. 

3. Provide realistic requirements for vibration isolation 
technology. 

The Space Shuttle Acceleration Environment 

The Orbiter Structural coordinate system shown in Fig. 

1. The accelerations experienced in a manned orbiting 
space laboratory are loosely classified as quasi-steady, 
oscillatory, and transient The quasi-steady classification, 
by convention, is assigned to accelerations with frequency 
below about 0.01 Hz. This encompasses accelerations due 
to atmospheric drag, gravity gradient effects, and rotational 
forces and is estimated [1] and measured [2,3] to be on the 
order of 1CT® g Q . The drag contribution varies with the 
orbital frequency of about 10"* Hz. The operation of 
machinery, control thrusters, solar panels, human activity, 
etc., contribute to higher frequency accelerations. 

Typical acceleration data taken on STS-47 (launched 
Sept 1992) are shown in Fig. 2 as a power spectral density 
plot This plot represents the environment during nominal 
crew activity and is characterized by distinct frequency 
peaks in the 1-10 Hz range that correspond to the Orbiter 
structural modes. The term nominal crew activity is used to 
indicate a time period when no significant acceleration 
sources such as vehicle maneuvers, water dumps, satellite 
launches and crew exercise are present The measured 
acceleration magnitudes are generally found to vary with 
the level of craw activity. Ftg. 3 shows an example of the 
difference in levels during periods of crew sleep, nominal 
crew activity, and crew exercise on a bicycle ergometer 
taken on the STS-40 mission [4], 

The data in each plot is a PSD calculation performed 
on 50 seconds of raw SAMS data and presented in units of 
Micro-g/Hz 1/2 . The RMS value for the data on each chart is 
shown in the figure. Acceleration levels in time vary from 
10* 4 to 10^g o . The SAMS sensor, located in rack 5 in the 
Spacelab module, recorded data at 25 sampies/sec with a 5 
Hz low pass filter. An important characteristic of the g-jitter 
environment is that the oscillatory and transient sources 
tend to excite Orbiter and payload structural modes. 
Because of this, the effects of a high magnitude transient 
source may be felt by an experiment in the form of 
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damped ringing for some time after the initial event 
occurred [4], The 3.5, 4.7 and 5.5 Hz components are all 
related to the excitation of structural modes by exercise 
activity and occasional thruster firings. Note the change of 
in the graphs. Even during the sleep period, the 
observation is dominated by the vehicle structural modes, 
but the general level has decreased by an order of 
magnitude. An example of the excitation of structural modes 
by an oscillatory source is the 17 Hz signal (not shown in 
Fig. 2). The KU-band communications antennae on the 
Orbiters dithers at -17 Hz to prevent stiction of the gimbal. 
This dither frequency tends to excite a 17 Hz Orbiter mode. 
During USMP-1, recorded not only the 17 Hz signal, but 
also first through fourth harmonics at 34, 51, 68, and 85 Hz. 
The most interesting feature of this signal is the extreme 
variability of its signature. Acceleration levels while the 
antenna was operating varied by an order of magnitude. 
The source of the variability is presently unknown, but it 
seems likely that it is related to the pointing angles of the 
system. 

Protein Crystal Growth 

Several proteins have been flown during past shuttle 
missions with the aim of growing bigger and better crystals. 
Protein crystal growth in microgravity benefits from the 
following advantages. 

1. The free suspension of the protein in the growth 
medium (due to diminished gravity) promotes uniform 
crystal growth and eliminates container interface effects. 

Z Convection effects are minimized and a diffusion 
controlled crystal growth can be realized. 

3. Internal stresses in the complex biological 
macromolecules are eliminated due to the reduced 
hydrostatic pressure environment This helps in improved 
internal order of the grown crystal and prevents the collapse 
of the big complex molecules. 

The data collected from past shuttle missions have all 
indicated that microgravity can be effectively utilized to 
obtain higher quality protein crystals. The space grown 
crystals tend to be significantly better than the best earth 
grown samples and thus are Ideal for protein 
crystallography. However, past experience has also shown 
that the microgravity environment on the Space Shuttle can 
sometimes yield unpredictable results in PCG. Poor crystal 
quality due to cracking, stunted crystal growth, multiple and 
some times showers of small crystals have been observed 
with certain PCG experiments. The causative factors for 
these effects have not all been completely understood. It is 
generally agreed however, that crew activity and g-jitter all 
play some role in this crystal degradation process. 

In any PCG experiment, the solution transport within 
the growth medium and the crystal surface attachment 
kinetics play key roles in determining the crystal growth 
rata. With most small molecule crystals, solution transport 
is the rate limiting step. In such situations, forced 
convection in the solution will increase the growth rate up to 
a certain point beyond which surface kinetics becomes the 


controlling factor. Experiments by Pusey et al. [6] using 
tetragonal lysozyme have shown that forced flow rates of 
30-40 pm/s slow and eventually stop the growth of 10 pm 
crystals. More recent experiments by Pusey [7] suggest that 
this growth cessation is present even at lower flow rates but 
the growth declines over a much larger time than at the 
higher flow rates. In the experiments [7], even very small 
velocities of the order of * 25 pm/s were seen to 
irreversibly terminate the (110) face growth of tetragonal 
lysozyme crystals. Corroborating results have been reported 
by Nyce and Rosenberger [8], who further showed that this 
growth cessation is not observed in the growth of 
orthorhombic lysozyme crystals. In their studies, Grant and 
Seville [9], proposed and subsequently disproved based on 
theoretical analyses a host of mechanisms that may 
account for the growth cessation behavior. The inclusion of 
contaminants from solution was hypothesized as the 
probable cause of this behavior. However, as mentioned in 
[7], this mechanism does not explain the absence of similar 
effects in the growth of orthorhombic lysozyme (8). Thus for 
some crystals at least, even small convective flows are 
deleterious for growth of the crystal. The exact mechanism 
for this effect is as yet not fully understood. 

Time Scales and Boundary Layers 

As a first step in the modeling process, estimates of the 
system time scale and boundary layer thickness have to be 
determined. Hens egg protein, lysozyme, is chosen for 
modeling purposes because it is a commonly available 
protein which is well characterized and and frequently used 
as a benchmark to test out new growth procedures, analysis 
etc. Estimates of its physical properties are shown in table 1. 


Table 1: Properties of Lysozyme 


Density, p, (g/cm 3 ) ^ 

1.2 

Density difference, Ap (g/cm 3 ) 

0.01 

Diffusivity, D (cnAs) 

icr 5 

Kinematic viscosity, v (cnAs) 

0.01 

Crystal Length, L (cm) 

0.1 

Crystal Radius, a (cm) 

0.05 

Grashof Number, Gr (1 g Q ) 

81.67 

Schmidt Number, Sc ( v/D) 

10 4 


The time scale of flow (fluid motion) can be estimated 
for both convective and diffusion dominated transport In 
terrestrial experiments, the flow is typically dominated by 
convection and the governing time scale is denoted as t c 
defined as 

x c = 8^/v (1) 

where 5 is the boundary layer thickness defined as 

5 = L(ScGr)* 1M (2) 

and the Grashof number is defined as 
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Gr = g(Ap/p)(L^/v2) (3) 

Therefore using the physical properties listed in table 1 
estimates can be drawn for 8 and x at different gravity 
levels. These are listed In table 2. 


Table 2: Time scales and Boundary layer extents 


g - level 

5 (cm) 

x (sec) 

1g 9 
1^3o 
10*? 9 0 
10- 3 g Q 

33.26x1 Or 4 

i.ixicr 5 

59.1 5x1 O' 4 

3. 50x1 0" 3 

105. 2x1 (T 4 

1.107x10- 2 

187.06x1 0 -4 

3.50x10-2 

icr 4 g„ 

0.03326 

0.1106 

icr 5 g Q 

0.0591 

0.350 


0.1052 

1.1079 


From table 2, it is evident that the boundary layer 
thickness is only of the order of 33 pm for earth grown 
crystals but is of the order of the crystal radius, which is 
incidentally also the diffusion length scale, at a g-level of 
10^g o . The length scale thus suggests that diffusion limited 
crystal growth can be achieved at a residual g-level of 10 p 
g. This remains to be verified by detailed fluid modeling. 
The diffusion time scale can then be evaluated as 

t- * a 2 /v = 0.25 sec. (4) 

w 

This implies that fluid flows can be initiated by excitations 
with time scales * 0.25 sec or by frequencies S 4 Hz. As 
di v- i Tcnri in an earlier section, the acceleration 
environment on the Space Shuttle contains several distinct 
vehicle structural modes in the 1-10 Hz range which can 
affect the fluid flows in protein crystal growth. The short 
time scale is due to the extremely small crystal size and as 
a result it takes but a small fluid velocity to traverse this 
dimension. One result of this fluid transport is the 
destruction of the existing solute gradient adjacent to the 
crystal. Thus an already established diffusional gradient is 
completely swept away by the fluid. The time scale for 
diffusion limited solute transport can be estimated as 

* a2/D = 2500 sec. (5) 

and thus it takes nearly 41 mins, for re-establishing the 
concentration gradient 

Order of Magnitude Analysis 

The magnitudes of the oscillatory disturbances that are 
of consequence to PCG can also be estimated by 
implementing a simple Order of Magnitude Analysis (OMA) 
of the equations describing the system. The analysis 
Involves choosing appropriate scaling factors for length, 
velocity, concentration, and other variables of importance to 
the experiment, while determining the dominant terms in the 
governing equations describing the system. Estimates of 
fluid sensitivity as a function of acceleration amplitude and 


frequency are obtained, and these estimates can be used 
as a preliminary guide to more detailed computations. A 
representative OMA curve for proteins is shown in Fig. 4 
1101 . 



FIG. 4 OMA CURVE FOR PROTEIN CRYSTAL GROWTH 


The analysis shows that the proposed space station 
microgravity environment would be conducive towards 
achieving this goal. PCG on board other manned vehicles 
on the other hand, would have to resort to other means to 
achieve the required acceleration environment Vibration 
isolation technology can be used to filter harmful 
frequencies and thus minimize the resulting convective 
effects. Current technology using active vibration devices 
can filter frequencies up to approximately 0.1 Hz and thus 
hold promise in this area. It should be pointed out that 
OMA estimates are based on a very simple model of the 
physical processes and probably over estimate the 
convective transport The model is however expected to 
provide correct trends and perhaps quantitative information 
within an order of magnitude accuracy. The level of 
accuracy can evaluated by detailed CFO modeling. 

Fluids Modeling 

The 2-0 computational model is shown in Fig. 5. A 
crystal of dimensions L = 2a » 0.1 cm is adopted for the 
calculations with the walls placed 5a and 4a to the side and 
below the crystal, respectively. Coding details are given in 
{5J. The crystal is modeled in planar two dimensions 
assuming an axis of symmetry. Comparisons of planar and 
axisymmstric calculations show that some differences exist 
in the two computed flow and solute fields. In an 
axisym metric situation, the effects of the walls are 
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attenuated whereas in a Cartesian system the wall effects 
are felt in some measure. The actual crystal growth process 
is typically of the order of 1-2 mm/day and hence is treated 
as a quasi-steady phenomenon. This implies that even 
though time accurate calculations are carried out to record 
the fluid response to time varying disturbances, the actual 
growth of a crystal phase is not modeled. This is a 
reasonable assumption considering the disparity between 
the respective time scales. No slip boundary conditions are 
applied on the crystal walls and on the computational 
domain x = L and y = 0. The understanding of the kinetics 
of the protein crystal growth process and its effect on 
surface flow characteristics is still in its infancy and hence 
the no slip velocity condition is imposed at the crystal wall. 

Residual-g Effects 

As a first step, the fluid response to terrestrial gravity 
(1g ) is calculated. An 81x51 (x, y) staggered grid is used in 
the calculation. Scaling analysis has shown that , the 
boundary layer thickness in this instance is about 33 pm. 
The grid spacing is chosen such that the distance from the 
crystal wall to the nearest grid point is 5.067 pm, so that 
there are about 6 grid points within the solute boundary 
layer. Results for this case are shown in Fig. 6. A 
maximum tangential fluid velocity of 70 pm/s is calculated 
adjacent to the crystal face. This compares well with an 
estimate of 90.35 pm/s obtained from scaling analysis 
where the maximum velocity is given by the expression 
(Gr/Sc) 1/2 (v/L). 

The equi-concentration contours or isomers for this 
case are also shown in Fig. 6. The steep concentration 
gradients next to the crystal are apparent from the figure. 


The boundary layer thickness is « 0.12 mm. The mass 
transfer into the crystal is characterized by the average 
Sherwood number defined as 

Sh = /(dC/an)*,,, (8) 

where n denotes the local normal to the wall. The Sherwood 
number is averaged over the entire crystal (L x 2a). For 1 
g 0 , the calculated Sh is 16.4, compared to Sh = 1.0 for pure 
diffusion limited growth. 

The fluid response to a parametric variation of the g- 
level (10*2g o , lO^g* 1<r*g 0 , lO^g,,. and 0g o ) is plotted in 
Fig. 7 as isomer distributions. The gradual thickening of the 
solute boundary layer is dearly visible as the g-level is 
decreased. The respective Sherwood numbers are also 
Indicated in the figure. The Sh values characterize the 
strength of convection within the cell and are seen to 
diminish as the gravity level is reduced. The diffusion 
limited case (g = 0) represents ideal crystal growth 
conditions and provides a reference point from which to 
compare the other growth scenarios. For example, at g - 
10*®g o (typical of the residual gravity level in space) only 
minor isomer distortions can be discerned from the figure 
and the Sherwood number is dose to unity. Diffusion 
limited growth conditions are approached at this low gravity 
level as predicted by scaling analysis. 

It is to be noted that the isomer contours for g * 0 show 
the effect of the computational domain by being drawn 
doser together near the bottom of the domain (below the 
crystal) than in the upper regions (above the crystal). Also, 
in order to maintain mass balance, the upward convection 
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plume in the system forces fluid entrainment and downward 
flow in the vicinity of the side wall. These effects are 
artifacts of the computational domain. The main interest 
however, is the fluid flow in the vicinity of the crystal itself. 
These axial flows are shown in Fig. 8 for the different 
residual g levels simulated. As the g-level is dropped, the 
velocity peaks are attenuated and shifted to the right as the 
boundary layer thickens. 
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Impulse type disturbances 

Transient events like thruster bums to achieve a certain 
flight attitude are more predictable than random Impulse 


type disturbances due to crew activity like opening and 
closing of lockers, the action of control relays to activate 
equipment or control experiments and astronauts bumping 
into equipment etc. Impulse strengths of up to ICT^q have 
been measured during past shuttle flights. In order to study 
fluids response to these disturbances, three specific events 
are simulated. 

Single Impulse (Ssturtance 

i. Ten seconds of zero gravity; Pure 0/L conditions. 

ii. One second impulse of strength; Case A: lO*^ and 
Case B: lO^g,,. 

ID. Calculations carried out for 600 seconds (10 minutes) 
after impulse with zero gravity. 

Fig. 9 shows the isomers and stream function plots at 
different times in the calculations for case A. The onset of 
the flow is immediate after the impulse and the flow 
development is complete 1 sec after the impulse. A 
clockwise rotational cell is estab lish ed in the system. No 
disturbance effects are however evident in the concentration 
field at this time. The maximum velocity right next to the 
crystal is about 4 pm/s and the maximum global velocity 
(entire computational domain) is about 40 pm/s. A time 
history of the global maximum velocity is shown in Fig. 10. 
There is a sharp spike in the flow which rapidly decays as 
shown in the figure. The transient 4 pm/s flow adjacent to 
as the crystal hardly disturbs the solute distribution because in 
(1 sec only a distance of 4 pm is traversed by the flow 
adjacent to the 100 pm crystal face. Some isomer 
disturbances are evident 100 secs after the impulse 
'"disturbance and dearly noticeable variations are seen 10 
mins after impact These variations are due to the 
drculation in the flow field set up by the impulse that lingers 
for a fairly long time due to fluid inertia and the absence of 
any significant damping effects. This drculation eventually 
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draws the isomers as shown in the figure. Fig. 10 also 
shows the maximum departure of concentration in the 
system from the reference case (g - 0) as a function of 
time. The disparity between the flow and solute responses 
is strikingly evident from this figure. In fact, for the solute 
Held, AC_« is still showing a growing trend 10 mins after 
the impulse. Results for an impulse of lO^g,, show similar 
results with the system response attenuated by an order of 
magnitude. 

Double Impulse disturbances 

The disturbance scenario for this part of the simulation 

is outlined below. 

L Ten seconds of zero gravity. Pure D/L conditions, 
ii. One second impulse of strength lO'^ggC+x direction), 
iil 1 second of zero gravity after impulse. 

iv. One second impulse of strength lO^g^x direction). 

v. Calculations carried out for 600 seconds (10 minutes) 
after second impulse with no gravity. 

A pulse/anti-pulse type disturbance is interesting 
because the mean g-level imparted to the system is zero 
and yet there exists net transport in the system response. 
Similar to the single impulse response, the flow responds 
almost instantaneously setting up a clockwise circulation 
after the first pulse and a complete reversal after the anti- 
pulse. Maximum velocities of 3.72 pm/s and 3.69 pm/s are 


recorded due to the first and second impulses respectively. 
Ten minutes into the computation, significant isomer 
distortions due to the residual flow circulation in the system 
is noticed. So, in essence, the flow has no memory of the 
first impulse but retains the memory of the second impulse. 
Similar results have been reported from computations in 
other configurations. 

OsaMory disturbances: 

The system response characteristics to an oscillatory 
disturbance frequency of 3.183 Hz (20 rad/s) and amplitude 
of lOr^o are shown in Fig. 11. This frequency is dose to a 
structural mode of 3.5 Hz. and also to the 2.5 Hz. noticed in 
SAMS data during treadmill operations on the Space 
Shuttle. The figures show complete flow reversal during 
each cycle. Fig. 11 shows the isomer and stream function 
plots during the 5th cyde of g-jitter and are representative 
of a periodic behavior attained within the system. The 
calculations show that the solute field is relatively 
unaffected by the disturbance but the flow field displays a 3 
Hz response as predicted by scaling arguments. This 
implies that a protein molecule is constantly being jiggled 
by the oscillatory disturbance. Whether such fluid motion is 
deleterious to crystal growth is undear. Time history plots 
of Sh and (AC) max also show frequency modulated 
response, but the variation magnitudes are very low, • 10* 5 . 
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Conclusions 


The salient results from the investigation are as follows: 

1 6-jitter dominates the spacecraft acceleration 
environment It is comprised of a myriad frequencies and 
displays no preferred orientation. The g-jitter magnitudes 

can be as high as 1 milU-g (10*^g o ). 

2. Impulse type disturbances are random in nature and 
hence unpredictable. The solutal field response to impulsive 
forces is especially long term and considerable. Impulse 
type disturbances are also deleterious to PCG in other 
respects e.g. drop dislodgment multiple crystals, crystal 
cracks etc. It is therefore prudent to take remedial measures 
to safeguard against their pernicious effects on matenals 

P 3° C pCG 9 observations and analyses indicate susceptibility 

^^Calculations show the PCG flow field is susceptible to 

1-10 Hz frequency range. . 

5 PCG is a good candidate for Vibration Isolation. 

Current vibration isolation system can filter up to 0.1 Hz. 
Passive systems can also filter down to this frequency level. 

8. A passive Isolation System for PCG Is a first 
recommendation. 
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